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宽带 大 规模 MIMO-OFDM 系统 自 适 应 稀 琉 信道 估计 方案 
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(太原 科技 大 学 电子 信息 工程 学 院 , 太原 030024) 


摘 要 : 大 规模 MIMO-OFDM 系统 下 行 链 路 利用 压缩 感知 算法 获得 信道 状态 信息 需要 已 知 信号 的 稀 玖 度 作 为 先 验 条 
件 ， 然 而 实际 环境 中 ， 无 线 信道 的 稀疏 度 是 未 知 的 。 利 用 大 规模 MIMO 信道 的 空 时 共同 稀疏 性 的 特点 与 不 同 SNR 下 
设置 不 同 停止 迭代 阔 值 的 思想 改进 压缩 感知 重 构 算 法 ， 目 的 在 于 使 所 提 算 法 不 仅 提升 估计 性 能 ， 还 可 以 准确 获得 信道 
的 动态 稀疏 度 。 通 过 实验 可 知 ， 相 比 传统 的 CoSaMP 算法 和 S-CoSaMP 算法 ，SSA-CoSaMP 算法 在 同等 信 骂 比 下 具有 
更 良好 的 信道 估计 性 能 ， 并 且 可 以 自 适应 地 获取 稀疏 度 。 更 适合 实际 工程 中 应 用 。 
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Adaptive sparse channel estimation for broadband massive MIMO-OFDM systems 


Huang Can, Li Suyuei, Wang Anhong 
(Collagel of Electronic Information Engineering, Taiyuan University of Science & Technology, TaiYuan 030024, China) 


Abstract: In massive MIMO-OFDM system, most of the compressed sensing theory reconstruction algorithms need the 
sparsity of the known signal as a priori condition. However, in the actual environment, the sparsity of the wireless channel is 
unknown. Therefore, the proposed algorithm used the space-time common sparsity of the Massive MIMO channel to 
reasonably set the stopping iteration parameters under different SNR to obtain accurate dynamic sparsity. The experimental 
results show that, compared with the traditional CoSaMP and S-CoSaMP algorithms, the SSA-CoSaMP algorithm has better 


channel estimation performance under the same SNR and can obtain the sparsity adaptively. 
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sparsity 
C o aë 家 现 阶段 的 蜂 帘 系 统 都 是 FDD 模式 的 加。 所 以 ， 对 于 FDD ft 


= 式 下 的 大 规模 MIMO 的 信道 估计 研究 具有 深远 的 意义 。 然 而 ， 
大 规模 MIMO 是 通过 利用 在 基站 侧 部 署 大 量 的 发 射 天 线 。 ”在 FDD 模式 下 ， 对 于 下 行 链 路 进行 准确 估计 是 十 分 困难 的 器 

且 实 现 同时 服务 于 多 个 单 接收 天 线 用 户 的 系统 中 。 因 此 ， 大 ”基站 侧 分 布 着 大 量 发 射 天 线 ， 而 这 些 发 射 天 线 同 时 向 用 户 发 送 
规模 MIMO 表现 出 展 好 的 抗 噪声 、 抗 窄带 衰落 的 能 力 ， 以 更 好 数据， 因此 ， 这 样 就 要 求 用 户 侧 准确 估计 出 每 一 根 发 射 天 线 与 
的 利用 日 渐 珍 贵 的 频谱 与 能 量 资源 ， 成 为 了 未 来 5G 发 展 的 关  ” 用户 侧 的 单 接收 天 线 所 构成 的 CSI。 同 时 满足 较 高 的 精确 性 和 


lm 


HERRA, (E AGNES MIMO 系统 中 ， 无 论 是 信号 检测 、 波 束 ，” 较 低 的 导 频 开销 也 是 十 分 困难 的 。 
赋 形 ， 还 是 资源 分 配 、 信 道 分 配 都 需要 获得 准确 的 CSI。 因 此 ， OFDM 系统 有 诸多 优势 ， 例 如 较 高 的 频谱 利用 率 ， 良 好 的 
性 能 杰出 的 信道 估计 方案 成 为 了 大 规模 MIMO 的 支撑 技术 。 抗 多 径 干扰 的 能 力 和 抗 窄带 衰落 的 能 力 等 由。 在 无 线 通 讯 用 户 


目前 , 关于 大 规模 MIMO 的 绝 大 多 数 的 研究 都 是 在 时 分 双 ”激增 的 今天 ， 需 要 提升 更 多 的 系统 容量 ， 因 此 把 正 交 频 分 复 用 
工 (time division duplex, TDD) 模 式 下 进行 的 , 目的 在 于 利用 上 下 ”系统 与 大 规模 MIMO 系统 联合 作用 , 从 而 实现 天 线 分 集 与 空间 
行 信道 的 互 易 性 ， 通 过 对 上 行 链 路 的 信道 估计 而 得 到 下 行 信道 复 用 中 。 当 前 的 通信 环境 下 ， 可 利用 的 无 线 带 宽频 率 越 来 越 医 
WI CSI, 但 是 ， 频 分 双 工 (frequency division duplex, FDD) 模 式 — Z5, 将 两 种 系统 结合 起 来 , 可 以 极 大 地 提升 频谱 资源 的 利用 率 ， 
比 TDD 模式 的 延 时 更 低 、 通 信 效 率 更 高 ， 且 国际 上 大 部 分 国 同时 还 可 以 针对 无 线 信 道 频率 有 选择 性 衰落 环境 进行 有 效 处 
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臣 型 的 宽带 无 线 传 输 环 境 存在 频率 选择 性 衰落 的 情况 ， 因 
此 ， 为 了 解决 这 一 问题 ， 我 们 将 大 规模 MIMO 与 OFDM 系统 
联合 使 用 建立 大 规模 MIMO-OFDM 系统 四。 


在 大 规模 MIMO-OFDM 传输 环境 中 ， 分 布 着 大 量 的 散射 
体 ， 形 成 了 瑞 利 多 径 传输 信道 。 并 且 无 线 信 道 的 绝 大 多 数 能 量 
关中 在 信道 冲击 响应 (channel impulse response，CIR) 的 极 少 分 
量 上 中 ， 使 得 用 户 侧 与 基站 侧 所 形成 的 CSI RRE T Fit 


性 [ 0,11] 


in 


o 


T 4 erst Pr 5 LER b PE, 可 以 利用 压缩 感知 (compressive 
sensing, CS) 对 信道 进行 估计 09。 然而 ， 实 际 环境 中 ， 无 线 信 道 
的 稀 玻 度 是 未 知 的 。 那 么 就 要 求 在 不 知道 稀 玻 度 的 情况 下 ， 算 
法 本 身 可 以 估计 出 信道 的 稀疏 度 031。 因 此 ， 本 文 提出 稀疏 度 自 
适应 压缩 采样 匹配 追踪 (SSA-CoSaMP) 算 法 , 利用 空间 共同 稀 玻 
性 不 仅 可 以 提升 估计 性 能 ， 还 可 以 节省 导 频 开销 ， 并 且 使 算法 
获得 稀 玻 度 自 适应 的 能 力 。 

本 文 的 主要 工作 包括 以 下 四 个 方面 

a) 对 CoSaMP 算法 进行 稀疏 度 自 适应 化 改进 ， 使 得 算法 在 
稀 玻 度 未 知 的 情况 下 可 以 自身 获取 稀疏 度 ; 

b) 在 稀 琉 度 自 适应 的 基础 上 利用 大 规模 MIMO 信道 空 时 共 
同 稀 琉 特 性 ， 进 行 结构 化 处 理 。 不 仅 使 获取 稀疏 度 速度 更 快 ， 
而 且 提 升 估计 性 能 及 准确 度 ; 

c) 算 法 停止 迭代 参数 a 的 确定 ; 

dj) 通过 仿真 程序 验证 了 提出 的 信道 估计 算法 的 六 


能 。 
1 大 规模 MIMO-OFDM RHE HAIE 


1.1 大 规模 MIMO 信道 的 空间 共同 稀 玻 性 
在 大 规模 MIMO-OFDM 系统 信号 传输 的 实际 环境 中 ， 分 
布 着 大 量 的 散射 体 ， 形 成 了 瑞 利 多 径 传 输 信道 。 并 且 无 线 信 
的 绝 大 多 数 能 量 集中 在 CIR 的 极 少 分 量 上 M4, 使 得 用 户 侧 与 
站 侧 所 形成 的 CSI CAU EL T PEU! 
与 此 同时 ， 基 站 侧 每 一 根 发 射 天 线 发 出 的 信号 在 到 达 用 
侧 的 接收 天 线 时 ， 所 经 过 的 散射 体 非常 相似 ， 几 乎 一 致 。 也 就 


dE uk 


Ir 


是 用 户 侧 的 单 天 线 与 基站 侧 的 所 有 天 线 有 空间 共同 稀疏 性 的 特 
点 。 此 外 ， 由 于 与 时 间 信 道 相关 的 路 径 增 益 相 比 ， 路 径 延 迟 的 


变化 慢 得 多 ， 所 以 在 相干 时 间 期 间 ， 这 种 稀疏 性 几乎 不 变 09， 
形成 时 间 相 关 性 的 特点 。 
内 型 的 单 小 区 下 行 多 用 户 大 规模 MIMO-OFDM 无 线 通信 


系统 存在 空间 共同 稀疏 性 ， 而 这 是 由 多 径 效 应 、 天 线 物理 特性 
和 无 线 传输 坏 境 等 多 方面 因素 共同 引起 的 。 在 无 线 传输 环境 中 ， 
散射 体 使 得 信号 从 同一 根 发 射 天 线 到 接收 天 线 之 间 形 成 多 条 路 
径 ， 形 成 多 径 效应 。 然 而 ， 就 在 这 些 传 输 路 径 中 ， 从 发 射 天 线 


到 达 接 收 天 线 的 时 间 也 是 不 一 致 的 "7。 因 此 ， 使 得 极 少 部 分 多 
径 分 量 上 有 着 CIR 的 绝 大 多 数 能 量 , 即 CIR 能 量 不 是 均匀 分 布 
在 信道 的 全 部 抽 头 上 ， 呈 现 稀疏 性 。 特别 地 ， 大 规模 MIMO 的 
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黄 W, 
第 m 个 发 射 天 线 和 一 个 单 接收 天 线 的 用 户 间 形成 的 CIR 可 以 表 
达 为 
hm [0,5 D) (D 


其 中 代表 第 r 个 OFDM. 符号 ，L 代表 等 效 信道 长 度 ， 
I1<m<M 。 
ITa ER h, 的 支撑 集 ， 可 表示 为 


Tp, =supp{h,,}= (&|h, [0] > aint ESL} 
其 中 : gs RRE P RE eME o RIE A L|, | 是 受 实际 
f TEI A e dE B PS A RE PE T A G, — A T UU 


m,r 


图 1 大 规模 MIMO 2s EH [o] fs or tk 


在 典型 大 规模 MIMO 中 ， 基 站 侧 的 天 线 阵 列 是 密集 型 的 ， 
天 线 间 的 间距 与 到 达 用 户 侧 单 天 线 的 信号 所 经 过 的 距离 相 比 是 


相当 小 的 ， 以 至 于 可 以 忽略 不 计 的 ， 并 且 每 一 对 收发 天 线 之 间 
的 信道 中 存在 公共 的 散射 体 。 因 此 ， 基 站 侧 的 每 一 根 天 线 与 用 


户 的 单 天 线 所 展现 出 高 度 相 似 的 路 径 延 迟 05， 如 图 1。 并 且 


不 同 发 送 天 线 对 之 间 的 CIR. 的 稀 玻 特性 绝 大 部 分 也 是 重合 的 。 


这 就 是 大 规模 MIMO 信道 的 空间 共同 稀疏 性 ， 即 对 于 第 r^ 


OFDM 符号 ， 不 同 发 射 天 线 保 持 相 同 的 稀疏 性 ， 可 以 表达 为 


D,,-I,,-..-Iy, (2) 


1.2 大 规模 MIMO 信道 的 时 间 相关 性 
在 时 间 方 面 ， 信 道 的 路 径 增益 变化 慢 很 多 Po， 


因而 显现 


O 
LI 


了 时 间 相 关 性 的 特点 ， 甚 至 是 在 快 时 变 的 情况 下 也 有 同样 的 特 


点 。 换 名 话说， 尽管 从 一 个 OFDM 符号 到 下 一 个 OFDM 符号 
路 径 增 益 变化 很 大 , 但 是 路 径 延 迟 在 连续 的 几 个 OFDM 符号 
几乎 保持 不 变 。 这 是 因为 时 变 信 道中 ， 


, 


H 


路 径 延 迟 的 变化 时 间 与 


系统 带宽 成 反比 ， 而 路 径 增 益 的 相干 时 间 与 系统 载波 频率 成 反 
比 上 。 也 就 是 说 ， 在 路 径 增益 的 相干 时 间 内 ， 因 为 路 径 延 迟 几 
乎 保持 不 变 ， 连 续 G 个 OFDM 符号 的 CIR 保持 不 变 的 共同 稀 


朴 性 ， 如 图 2 Hrs. B 


REPHSERESIM 


相干 时 间 内 连续 的 不 同 OFDM 符 号 


空 时 共同 稀 朴 性 示意 图 


m,r+G-1 (3) 
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录用 稿 


综 上 ， 将 空间 共同 稀疏 性 和 时 间 相 关 性 合 称 为 时 空 共同 稀 


RELA, 


2 大 规模 MIMO-OFDM 系统 模型 


建立 一 个 典型 的 FDD 模式 大 规模 MIMO-OFDM 系统 ， 小 
区 内 基站 侧 天 线 数 为 M， 单 天 线 用 户 数 为 K， 且 M EK, "nl 
3 所 示 。 其 中 ，OFDM 符号 的 子 载波 总 数 为 N， 基 站 第 m 根 天 
线 发 射 的 信号 (包含 数据 和 导 频 ) 为 S, e C" (m-12,., M) « H 


么 ， 用 户 侧 接收 到 的 信号 可 表达 为 
y- > S,Fh,+n (4) 
CP i. 
[oo] y [n.o 
terre "M 
zl. -1 J= ER Ec [yi i — i] 
. CP . 
。| | | 并 M 
at [no] y m r*a], 
; EE E: 
adn ki] FA UR 3 m Yn & - 1]. 
图 3” 大 规模 MIMO-OFDM 框架 
除去 保护 间隔 和 进行 DFT 变换 后 , 用 户 侧 根据 导 频 分 配方 
案 s 从 接收 到 的 信息 y 中 提取 出 第 r 个 OFDM 符号 中 的 导 频 信 


© y, , y, 可 以 表达 为 


y,- DX VFI; 


h, 
ln i +N, Ey C,F, |, ^, +n, 
On ya 


m-l 


- Soh, „+n, (5) 


m-l 


其 中 : 5 表示 在 一 个 OFDM 符号 中 ， 导 频 所 占 子 载波 的 位 置信 


Bm 


本 文 导 频 设计 方案 为 非 正 交 化 导 频 分 配方 案 ( 详 见 第 3 
T)O,-C,F|,: C,-diag(o,]: c,eC*" Xm m HEX 
线 的 导 频 幅 值 的 集合 。N, 表示 在 一 个 OFDM 符号 里 的 导 频 数 ， 
即 N, =ó FERKA NxN 的 DFT 4ER, F 是 取 握 的 
LIT F|, 表示 根据 8 BUE, TIER, n 表示 第 r 个 OFDM 
符号 在 信道 中 的 加 性 高 斯 白 噪声 向 量 。 

将 式 (5) 化 简 可 得 


y,= Dh +n, (6) 


其 中 DB=[D,D,,..., Dy] eC 


h, Es [rale m,] € pu B 


在 大 规模 MIMO-OFDM 系统 中 ，N, < ML。 因 此 式 (6) 是 
个 欠 定 方程 ， 用 传统 的 信道 估计 方法 ， 如 最 小 二 乘 (least 
square, LS) 算 法 不 能 可 靠 的 进行 估计 。 但 是 ， 经 过 上 述 的 分 析 


可 知 户 是 一 个 高 维 稀 玻 信 号， 因此， 可 以 利用 压缩 感知 理论 从 
HERO y, 将 所 重 构 出 来 。 
此 外 ， 利 用 1.1 节 所 提 到 的 大 规模 MIMO 固有 的 空间 共同 
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稀 玻 性 可 以 提高 CS 算法 的 重 构 性 能 。 


4 d, -[h, (0... n, (O)] ， 进 而 对 所 进行 重新 排列 ， 可 


以 得 到 : 
9,=[ 9,,0,,...,9!,| ecu (7) 
Ho 理 ， ^ o P[O, pL, Dp] 
[is 和 sw]sCww%， 那 么 由 进行 重新 排列 后 ， 可 以 表达 
为 
V =[P, PnP] E Cn (8) 
Kb, IR(8) ET DERE 
y, - d, +n, (9) 
经 过 上 述 利用 无 线 大 规模 MIMO ALIIS EREET 
的 转换 变形 , g 的 稀疏 性 得 到 结构 化 处 理 , 使 稀疏 性 得 到 增强 。 


利用 12 节 提 到 的 时 间 相 关 性 ， 大 规模 MIMO 系统 的 空间 
k 同 稀疏 性 在 相干 时 间 的 几 个 连续 的 OFDM 符号 内 是 几乎 不 
变 的 。 因 此 ， 在 导 频 分 配方 案 S. 下 ， 对 相干 时 间 内 连续 的 R 个 
OFDM 符号 进行 处 理 ， 式 (9) 可 以 整理 为 
Y-VH«n (10) 
Y=[ AY av agalet t 


其 中 ， 


Hagen Hs dd ics d, 1| eC a 


H gg JRRI ERAS T H-[HLHL...HTT 3 


利用 无 线 大 规模 MIMO 信道 的 时 间 相 关 性 ， 在 (9) 的 基础 
上 ， 进 一 步 使 信道 CIR 矩阵 结构 化 得 到 (10)， 提 升 稀疏 性 ， 更 
好 的 提升 CS 算法 的 估计 性 能 


3 ” 非 正 交 化 导 频 设计 方案 


导 频 对 于 信道 估计 有 着 至 关 重 要 的 作用 ， 优 秀 的 导 频 设计 
方案 可 以 极 大 地 提升 信道 估计 的 性 能 ， 如 果 设 计 的 导 频 方案 不 
合理 或 者 不 适合 当前 系统 ， 不 仅 不 能 进行 准确 的 信道 估计 ， 
至 增加 导 频 开销 ， 浪 费 宽 带 资 源 和 能 量 。 
31 正 交 化 与 非 正 交 化 导 频 

传统 正 交 化 导 频 设计 方案 是 在 奈 奎 斯 特 采样 定理 的 框架 下 
设计 的 ， 它 要 求 不 同 天 线 上 的 导 频 占据 不 同 的 子 载波 ， 并 且 天 
线 上 已 安排 了 导 频 的 子 载波 在 其 他 天 线 相 对 应 的 位 置 上 被 要 求 
安排 成 空 导 频 。 现 下 ， 正 交 化 的 导 频 设计 方案 已 经 被 广泛 的 使 
用 在 现存 的 MIMO 系统 中 ,因为 天 线 数量 很 少 ， 所 以 导 频 开销 
不 是 很 严重 的 问题 ， 如 8 天 线 的 LTE-advanced 系统 。 但是， 当 
发 展 到 大 规模 MMO 时 , 基站 侧 天 线 配置 量 达到 128 根 乃 至 更 
多 时 ， 这 个 问题 就 不 能 再 忽视 了 ， 导 频 开 销 会 非常 严重 ， 造 成 
极 大 的 能 源 和 频谱 浪费 。 
基于 上 述 问题 ， 我 们 根据 文献 [23]， 在 本 文 应 用 了 非 正 交 
化 导 频 设计 方案 , 它 是 基于 CS 理论 进行 设计 的 。 与 传统 正 交 化 
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导 频 方案 不 同 ， 它 允许 不 同 天 线 上 的 导 频 完全 占据 相同 的 子 载 


波 。 通 过 利用 大 规模 MIMO RAA RE, HTE 


计 的 导 频 开销 可 以 实质 性 的 降低 。 
5 1 al: B ns 
* * »- 无 效 空 
|^ BSBH 
: po po m mea 
in nose 
ji i TE CEN 
Tx! Tx2 Tx Uma TM RS Tx1 Tx2 Tx M RŠ 
图 4 正 交 化 和 非 正 交 化 导 频 设计 方案 


3.2 非 正 交 化 导 频 设计 理论 推导 


在 CS 理论 中 ,根据 文献 [23]， 式 (10) 中 感知 矩阵 风 的 设计 


对 于 可 以 高 效 可 靠 压缩 高 维 稀疏 信号 是 非常 重要 的 。 在 大 规 
模 MIMO 信道 估计 中 ，w 只 由 导 频 分 配方 案 3 和 导 频 序列 
{G6,}" 决定 确定 ， 所 以 对 由 的 设计 可 以 转化 成 为 对 5 和 
(c,) 的 设计 。 根据 CS 理论 可 知 ， 如 果 风 具有 列 之 间 的 相关 


性 ， 对 于 进行 可 靠 的 恢复 稀 玻 信号 是 非常 有 帮助 的 。 这 也 就 局 


发 了 对 6 和 {6 oh 的 设计 思路 。 


首先 从 导 频 序列 {c y 设计 的 角度 着 手 研究 , 目的 在 于 使 


18 V, 中 的 任意 p 列 之 间 具 有 的 互相 关 性 


只 由 {c,j 本 身 决 定 ， 即 


VA F ) V, zu »y v 


-(s, «F^. ^) 


其 中 : F =F]; > I1<m<m SM 。 


较 小 。 而 这 种 互相 关 性 


-(c, c (1) 


为 了 获得 较 小 的 (yw we, 


, 


> 


DRI (0, Ve" 遵循 


=], 


独立 同 分 布 的 均匀 分 布 U[0,2x) ， 其 中 giten 表示 c ec 的 


第 k 个 元 素 。 并 且 w 的 每 一 列 的 


Lf, 范 数 都 为 常量 ， 即 


V uc 2 = JN, iz 此 时 可 以 得 
Kw. E V, (c ) C 
lim 1 as* | _ im T m 0 (12) 
Np >00 A Vp, t 5 N, >% N, 


根据 随机 和 矩阵 理论 (random matrix theory, RMT), HÆ 
实际 状况 下 的 N, 是 有 限 的 ， 因 此 根据 式 (12) 可 得 ， 非 正 交 化 导 


p 


频 序列 (c M. 可 以 使 得 U, 中 的 任意 
关 性 。 


列 之 间 具 有 良好 的 互相 


其 次 从 导 频 分 本 方案 5 设计 的 角度 继续 深入 研究 wy, ,与 


y,,, ,, 的 互相 关 性 。 在 0 Le 的 情况 下 ， 


Vt x 


小 的 |(w 


Bt 
E 


型 的 大 规模 MIMO 系统 中 ， 


如 ， 在 典型 的 基站 侧 配 置 128 根 天 线 的 大 规模 MIMO 系统 中 


6 的 设计 也 可 以 获得 较 


FEN, > 工 的 特点 。 例 


每 根 天 线 上 至 少 要 有 一 个 导 频 才 可 以 保证 该 天 线 可 以 进行 信道 
"i. D GNE QUNM DN 此 外 ， 信 
道 最 大 的 时 延 拓 展 是 3 ~ 5ys 并 且 典 型 的 LITE-Advanced 系统 的 
带宽 是 10MHz， 因 此 p «128 ^. EIE, SIA N, 必须 大 
于 工 。 
基于 N, > 工 这 个 条 件 , 如 果 用 均匀 间隔 的 方式 进行 导 频 分 


> 


c 
dr 


配 不 仅 使 用 范围 更 宽泛 ， 而 且 可 Lb 


f8 (v, ,) v, | Sv. 


匀 导 频 间 隔 为 


六 | mw 有 we 的 内 入 


L P 


Vous (e) qe 


slas F9) (e, « Ff?) 


N H 


: . Am : 
= Pep (jE + jÓ, ) 


k=1 
v .2m : 
x> exp| j —— £1 (Kk) + jO, n 
per! N Ma 
Yew[ i^ i09» int, "ES 


Koh: gy-r 是 导 频 子 载波 的 索引 集 ， 


1<l=l, -4 SL-1 > 0,70, -Oem © IE (rp 是 根 


SEEELLLIPA 


p 


| ARA paan) 中 得 到 的 ， 


N H se ; 
I(k) 5 I, +(k (x 1<k<N,， 其 中 : 万 是 第 一 个 导 频 
p 


的 子 载波 集合 且 1<7 «x 


因此 ， 式 (13) 也 可 以 表达 为 


XE E i. -(k E [eaan] (14) 
GVA N, 


term Ce e[0,1) ), 进一步 对 式 (14) 进 行 整理 可 


(Wn A y V, t, 


N 
p 2 p 
=h ev [2 Z -e)na (15) 


p 


考虑 当 m =m 时 的 情况 ， 此 时 有 A6，, =0， 则 式 (15) 可 化 
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H < 2z- 
(v... v. ner ( 5 0 - ne) (16) 
k=l 
Job. qo 人 <1 表 示 导 频 占 比 ， 因 此 当 ne =0 ， 此 时 ， 
可 以 得 到 
H 
lim (v...) Vt 
N, >% 
p 
r. ] — e1270 ne) 
= lim ii - NP 
p j£ t-e) 17 
[^ "E | Ha 
22m 2-5 eme 
其 中 : 村 1<7<L-1 和 L«N,; 则 ZA IS NLIS 


了 式 (17) 的 合理 有 效 性 。 
考虑 当 m zm 时 的 情况 ， 则 式 (15) 可 以 表示 为 


N, m 
(VAR M Wan 2 exp( Jő ) (18) 


k,m 


Hh. A-U 遵循 独立 同 分 布 的 均匀 分 布 


U[0,22) H1«k « N, 。 类 似 于 式 (12)， 可 以 得 到 


lim Uma) Vs 
N, >% 
P 
N, z 
i0, 
= lim Low ) =0 dy 
N,—00 


p 


综 上 所 述 ，N, 有限 的 实际 情况 下 ,根据 RMT， 对 于 wy 的 
任意 两 列 间 ， 式 (12)(17)(19) 的 渐进 正 交 性 表明 了 非 正 交 化 导 频 
分 配方 案 5 和 导 频 序列 (o Vt 都 获得 了 很 好 的 互相 关 性 。 此 
外 ， 在 基于 CS 理论 的 信道 估计 技术 中 ， 相 较 于 广泛 使 用 的 传 
统 随 机 导 频 分 布 方案 2， 均 匀 间 隔 导 频 分 布 方案 在 现实 环境 
更 容易 实施 。 更 重要 的 是 ， 均 匀 间 隔 导 频 分 布 方案 已 广泛 应 用 
在 现 阶段 的 蜂 帘 网 络 中 , 因此 它 使 大 规模 MIMO 系统 可 以 更 好 
地 向 后 兼容 进 现 阶段 蜂 帘 网 络 。 


4 ”改进 方案 


41 稀 琉 度 自 适应 原理 
在 实际 通信 环境 中 ， 信 道 的 稀疏 度 是 未 知 的 ， 并 且 信 道 状 
态 随时 间或 者 空间 的 变化 而 变化 。 因 此 ， 为 了 更 好 的 可 以 在 工 


程 中 应 用 ， 我 们 以 阀 值 思想 为 基础 ， 提 出 动态 稀疏 度 自 适 应 改 
进 方案 。 不 同 于 传统 重 构 算法 根据 稀 疏 度 确定 迭代 次 数 ， 该 方 


案 是 通过 设 定 阔 值 来 确定 迭代 次 数 的 ， 即 当 残 差 和 观测 值 之 间 
处 于 停止 欠 代 参数 所 要 求 的 比值 关系 时 算法 停止 迭代 。 并 且 ， 
停止 欠 代 参数 并 不 是 定 值 ， 它 在 不 同 信道 状态 下 ， 停 止 迭 代 参 
数 也 不 同 。 信 道 状态 我 们 用 SNR 来 表征 。 据 此 ， 可 表示 为 

IV <2 xl (20) 
Hob: V 表示 残 差 ，Y 表示 观测 值 。 表示 停止 迭代 参数 ， 它 
不 是 一 个 固定 的 值 ， 它 在 不 同 SNR 下 是 不 同 的 。 
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REDERI SUARE, MATERE 
KERTAA ETE 77; OF LEE 8 p 的 动态 设 定 ; 不 同 的 信 
道 状 态 下 ， 噪 声 对 观测 值 的 影响 也 会 有 不 同 ， 因 此 停止 欠 代 参 
数 p 不 能 设 定 为 一 个 固定 的 值 ， 而 应 该 在 不 同 的 信道 状态 下 是 
不 同 的 ， 用 SNR 表征 信道 状态 ， 也 就 是 说 ， 不 同 SNR 下 ， 售 
止 和 迭代 参数 8 是 不 同 的 。 
在 文献 [26] 中 也 提 及 了 利用 阔 值 的 方法 来 使 算法 可 以 自 适 
应 获取 稀 疏 度 ， 但 设 定 的 停止 迭代 参数 是 定 值 ， 不 能 根据 信道 
状态 变化 停止 选 代 参 数 ， 而 本 文 提 出 的 算法 可 以 动态 获取 稀 玻 


度 。 
4.2 考虑 稀疏 度 自 适应 的 大 规模 MMO 压缩 感知 信道 估计 方 
案 


利用 稀 玻 信道 空 时 共同 稀疏 性 和 稀 玻 度 在 未 知情 况 下 可 上 
身 获取 的 特点 ， 从 结构 化 和 动态 稀疏 度 自 适应 两 个 方面 提升 其 
EERE, PEETRE BEM CoSaMP (SSA-CoSaMP) 算 法 。 
表 1 SSA-CoSaMP 算法 流程 


SSA-CoSaMP 算法 具体 如 下 : 


输入 : 信号 观测 值 yo MWEE gy ， 最 大 信道 长 度 L， 基 站 侧 配 置 天 线 数 


M， 连 续 OFDM 符号 数 R， 算 法 停止 参数 B 。 


输出 : 估计 的 CIR HERE f o 


a) 初 始 化， 初始 残 差 V =Y; WAE S=1:; 循环 次 数 大 =1 ; LER 
Q, =% 


bz-uy 


c) 利 用 大 规模 MMO 1825 [| fibi PEE Gt P RE REGE epe, BU 


c(t) m POM 


c(z) -[c(0),....c(L-1)]| > zi" zH mtr nl. 


2 


0<zr< 工 -1 Q 其 中 


9 Q-QUsupp(o),,] 合并 集合 


> 


9r-OoU[O- L]U....U[O.-- L(M —D)] 合 并 集 


DV, = y — P” A 对 余 量 进行 更 新 ; 


DEWE |V] < DBxlY WES EWS k=k+l 8 =S+1 转 步 


Jk b). 


从 改进 后 的 SSA-CoSaMP 算法 流程 可 知 , SSA-CoSaMP 算 
法 并 不 是 根据 提前 给 定 的 稀疏 度 来 停止 迭代 ， 而 是 建立 残 差 和 
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观测 矩阵 间 的 关系 ， 当 满足 设 定 条 件 时 ， 自 动 停止 迭代 ， 从 而 
实现 了 算法 可 以 自 适 应 获取 稀疏 度 ， 为 该 算法 在 实际 的 信道 估 
计 中 的 工程 应 用 提供 可 能 。 从 SSA-CoSaMP 流程 可 知 ， 主 要 分 
为 两 部 分 ， 第 一 部 分 是 对 稀 玻 度 s 的 确定 。 算 法 在 初始 化 阶段 
WE ME s =1, 然后 利用 循环 友 代 的 方法 最 终 得 到 满足 停止 
ARERR s> ZI ERRAR aM h). BRE 
EMERE s 的 情况 下 , 利用 无 线 大 规模 MIMO 信道 的 空 时 
# 同 稀 朴 性 特点 ， 在 对 p 进行 重 构 估计 的 过 程 中 ， 每 次 迭代 都 
可 以 同时 对 多 个 稀疏 向 量 进行 更 新 , 对 应 于 SSA-CoSaMP 算法 
中 的 步骤 b)~g)。 最 终 得 到 准确 的 稀 纹 度 和 信道 CIR. 

如 果 对 CoSaMP 算法 只 利用 空 时 共同 稀疏 性 而 不 考虑 稀 玻 
度 自 适应 改进 的 话 ， 我 们 将 其 称 为 结构 化 CoSaMP (HJ 
S-CoSaMP) 算 法 。S-CoSaMP 算法 只 能 根据 给 定 的 稀疏 度 进行 
循 坏 迭 代 的 次 数 ， 而 不 是 由 实际 的 信道 状态 进行 判断 是 否 结束 
算法 。 这 样 带 来 的 后 果 有 两 种 :a) 算 法 已 经 重 构 出 结果 ， 但 算法 
本 身 为 了 满足 所 给 定 的 稀 玉 度 却 还 需 继续 循环 迭代 ， 造 成 时 间 


m 
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为 非 正 交 化 导 频 分 配方 案 ， 噪 声 由 零 均 值 单位 方差 的 复 高 斯 随 
机 变量 组 成 。 系 统 模型 为 国际 电信 联盟 车 载 BITU-B) 信 道 模 
70, 仿真 中 , 还 对 所 提出 的 SSA-CoSaMP、S-CoSaMP 、CoSaMP 
算法 进行 了 比较 。 在 性 能 的 衡量 方面 ， 用 归 一 化 均 方 误差 
(NMSE) 来 作为 指标 .NMSE 越 小 , 估计 出 来 的 结果 误差 就 越 小 ， 
则 表明 估计 性 能 就 越 好 。 

在 本 章 中 , 均 方 误差 归 一 化 表示 ， 即 NMSE, 并 如 下 定义 ， 
来 作为 信道 估计 性 能 的 衡量 指标 。 


— 


--l gl&-nf 
NMSE= E| -H| Q1 


作为 对 比 的 算法 ，CoSaMP 算法 参考 文献 [27]，S-CoSaMP 
算法 参考 文献 [28] 。 
5.1 停止 迭代 参数 5 的 确定 
随 着 信 品 比 的 不 断 降低 ,噪声 对 信号 观测 值 Y 的 影响 将 会 
增 大 ， 因 此 在 不 同 信 品 比 下 ， 也 会 相应 变化 。 本 节 从 对 比 思 
想 的 角度 考虑 ， 利 用 Exact LS 算法 来 获取 p 在 不 同 SNR. 下 的 


和 能 量 的 浪费 ，b) 算 法 重 构 出 的 结果 还 不 够 精确 ， 但 是 给 定 的 
稀 玻 度 已 经 满足 循环 途 代 的 次 数 ， 和 迭代 停止 ， 造 成 估计 的 精度 
不 足 。 在 接 下 来 的 第 6 节 仿 真 数据 及 结果 分 析 中 ， 给 出 了 
S-CoSaMP 和 SSA-CoSaMP 算法 的 性 能 对 比 。 

SSA-CoSaMP 算法 中 g 在 不 同 SNR. 下 是 不 同 的 .如 果 p 的 
值 是 固定 的 ， 而 不 是 根据 信道 状态 而 动态 变化 的 会 直接 对 重 构 


结果 造成 极其 严重 的 影响 。 若 9 设置 的 过 小 ， 那 么 会 使 算法 的 
迭代 次 数 增多 ， 进 而 引起 算法 用 时 过 长 ， 不 利于 工程 应 用 。 相 
反 ， 如 果 p 设 值 的 过 大 ， 那 么 会 造成 算法 帮 代 过 早 结束 ， 重 构 
出 的 结果 误差 过 大 ， 降 低 了 算法 的 重 构 精度 。 所 以 ， 的 设置 
对 重 构 精 度 和 迭代 时 间 都 有 至 关 重 要 的 影响 。 因 此 ， 在 接 下 来 


的 第 6 节 中 ， 我 们 通过 仿真 实验 确定 参数 p, 
精度 的 同时 将 可 能 的 降低 运算 时 间 。 


5 ”仿真 结果 及 分 析 


大 规模 MIMO 无 线 信 道 估计 技术 主要 以 均 方 误差 (MSE) 来 
评价 其 性 能 。 可 以 从 以 下 三 个 方面 试验 来 体现 : a) 估计 精度 对 
比 实验 ,估计 出 的 多 径 分 量 与 实际 建 模 相对 比 , 用 直观 的 方法 来 
表征 估计 性 能 :b) 随 导 频数 变化 的 性 能 仿真 实验 ， 导 频数 的 降 
氏 必 然 会 影响 估计 性 能 ， 用 导 频 数 变化 来 表征 估计 性 能 可 以 更 
好 的 反应 估计 算法 是 否 能 在 降低 导 频 开销 的 同时 又 不 过 多 的 影 
响 估 计 性 能 ，c) Bü SNR 变化 的 性 能 仿真 实验 ，SNR 可 以 很 好 
也 反映 信道 状态 ,用 SNR 变化 来 表征 估计 性 能 可 以 很 好 地 反映 
出 估计 算法 在 不 同 信道 状态 下 的 估计 情况 。 所 以 本 方案 主要 从 
以 上 三 个 方面 的 实验 来 说 明 方案 的 优越 性 。 

本 节 利 用 MATLAB 仿真 软件 对 所 提出 的 大 规模 
MIMO-OFDM 无 线 通 信 系 统 中 的 SSA-CoSaMP 算法 进行 性 能 
验证 。 系 统 参数 设 为 基站 侧 发 射 天 线 数量 为 128 根 ， 系 统 中 系 
统 子 载波 数 为 1024， 信 道 多 径 数 为 256， 连 续 OFDM 符号 数 
R=5 且 每 个 信道 为 独立 同 分 布 的 瑞 利 衰落 信道 ， 导 频 分 配方 案 


使 之 在 满足 重 构 


a 


1B 

言 道 估计 算法 一 般 可 分 为 两 个 阶段 ， 第 一 阶段 ， 估 计 重 构 
出 抽 头 所 在 的 位 置 ;第 二 阶段 ,估算 立 置 抽 头 的 值 。 而 Exact 
LS 算法 直接 给 定 了 抽 头 所 在 的 准确 位 置 , 只 需 估 算出 该 位 置 抽 
头 的 值 即 可 。 故 而 Exact LS 算法 代表 了 CS 类 估计 算法 的 性 能 
基准 。 因 此 ， 在 本 节 实 验 中 , 将 8 逐渐 变化 的 SSA-CoSaMP 算 
法 与 Exact LS 算法 进行 NMSE 对 比 。 当 恰好 使 两 者 NMSE 相 
等 的 8 ， 就 是 该 SNR 下 的 取 值 。 


k 


Em 
xS 


EE 
x 


NMSE 


(a)SNR-10dB 


一 白 ExactLS 算 法 
一 一 Xx SSA-CoSaMP t.i? 


(b) SNR=20dB 
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ve 
录用 稿 黄 h, $: 宽带 大 规模 MIMO- os 系统 自 适 aX ERI y 
A | 
停 正 参数 B i 
(c) SNR=30dB 
图 5 不 同 SNR FSER B (b) N, =900 
图 5(a)rP, g =0.058 时 ， 两 者 的 NMSE 恰好 相等 。 当 大 于 图 6 不 同 的 信道 估计 算法 的 准确 性 对 比 (SNR=20dB) 
0.058 时 ,两 者 NMSE 不 一 致 。 当 小 于 0.058 时 , 虽然 两 者 NMSE — 5.3 信道 估计 算法 随 导 频数 变化 的 性 能 仿真 
相等 ， 但 不 是 恰好 达到 相等 的 临界 取 值 。 导 频 是 在 基站 侧 被 安排 在 OFDM 符号 某 部 分 特定 子 载 ; 
同 理 , 根据 图 Sb) 和 图 5(c) 的 实验 仿真 结果 , 可 以 归纳 出 8 上 的 特殊 数据 ， 接 收 端 根据 收 到 的 导 频 来 对 信道 进行 估计 。 导 
取 值 范围 ， 如 表 2 所 示 。 频数 据 占用 了 一 定 的 传输 资源 。 为 了 提高 系统 的 频谱 效率 ， 必 
表 2 与 SNR 相对 应 的 算法 停止 参数 然 要 降低 导 频 数据 在 每 一 帧 OFDM 符号 中 的 比例 , 而 降低 导 频 
SNR 对 应 的 g 数 势 必 会 影响 估计 性 能 。 因 此 ， 进 行 导 频数 变化 对 信道 估计 算 
<10dB 0.058( 图 5(a)) 法 性 能 影响 的 仿真 试验 是 十 分 必要 的 。 
10dB 一 25dB 0.054( 图 5(b)) ERT 
>25dB 0.052( 图 5(c)) x CIERRE 
5.2 信道 估计 算法 的 精度 对 比 实验 中 


对 算法 估计 出 结果 的 优 劣 判断 可 以 从 两 个 方面 来 考察 ， 即 
十 计 出 抽 头 的 位 置 和 幅 值 。 因 此 ， 本 节 将 算法 估计 出 的 结果 与 
实验 所 构建 的 信道 模型 进行 以 上 两 个 方面 的 对 比 来 验证 其 性 


fE o 
图 6(a) 是 在 SNR=20dB 且 导 频数 为 700 时 进行 试验 的 仿真 EN | | | 
图 。 根 据 实 验 结果 可 知 ，CoSaMP 算法 估计 出 的 结果 无 论 是 抽 pe 
头 的 位 置 还 是 幅 值 都 不 准确 ， 而 S-CoSaMP 算法 的 结果 也 只 图 7 提出 算法 随 导 频数 变化 的 NMSE 性 能 对 比 
部 分 抽 头 的 位 置 准确 ， 其 他 结果 的 抽 头 位 置 和 幅 值 都 不 准确 。 7 是 SNR 为 10 dB 和 20 dB 时 ，SSA-CoSaMP 算法 、 
相 比 之 下 ，SSA-CoSaMP 算法 准确 了 很 多 ， 估 计 出 的 抽 头 无 论 。 S-CoSaMP 算法 和 CoSaMP 算法 的 NMSE 随 导 频数 增加 而 变化 
位 置 还 是 幅 值 都 很 准确 。 的 仿真 图 。 从 实验 结果 可 以 得 出 ,在 不 同 SNR 下 ， 随 着 导 频 数 
当 SNR=20dB ， 导 频数 为 900 时 ， 我 们 将 算法 估计 出 的 结 ”的 不 断 增 加 ， 三 种 算法 的 估计 性 能 都 有 所 改善 。 当 导 频 数 低 于 


Ke 

果 和 实验 建 模 进 行 了 对 比 仿真 ， 可 以 得 到 图 6(b)。 根 据 实验 结 ”750 时 ，SSA-CoSaMP 算法 的 性 能 最 优良 ， 其 次 是 S-CoSaMP 
果 可 得 出 ， 此 时 CoSaMP 算法 仍然 不 能 进行 准确 的 估计 ， 而 ”算法 , 而 CoSaMP 算法 的 估计 性 能 最 差 。 当 导 频 数 高 于 750 时 ， 
SSA-CoSaMP 算法 和 S-CoSaMP 算法 的 估计 性 能 基本 一 致 ， 都 S-CoSaMP 和 SSA-CoSaMP 算法 的 估计 性 能 基本 一 致 。 
非常 准确 。 5.4 ”信道 估计 算法 随 SNR 变化 的 性 能 仿真 

除了 第 5.3 节 仿 真 的 导 频 数 变 化 对 估计 算法 的 性 能 有 所 影 
响 外 ， 还 有 一 个 参数 对 估计 算法 的 性 能 也 有 很 大 的 影响 ， 那 就 
是 SNR。 为 了 进一步 验证 算法 的 性 能 与 SNR 的 变化 之 间 的 关 
系 , 所 以 本 节 进 行 了 信道 估计 算法 随 SNR 变化 的 性 能 仿真 的 实 
验 。 图 8 是 信道 估计 算法 随 SNR 变化 的 NMSE 仿真 图 。 并且 
分 别 选 取 了 导 频 数 为 650、750、850 时 的 情况 ， 使 实验 结果 更 
RE | 清晰 准确 。 

ms. 从 图 8(a) 可 以 得 知 ， 当 导 频 数 为 650 时 ， SSA-CoSaMP 算 
(a) N, =700 法 的 性 能 比 S-CoSaMP 算法 优秀 了 许多 。 而 CoSaMP 算法 随 着 
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SNR 的 增加 并 没有 明显 变化 ， 处 于 无 法 进行 出 有 效 运算 的 状 


态 。 该 实验 结果 表明 了 改进 后 的 算法 在 低 导 频 的 情况 下 ， 仍 能 
很 好 的 估计 出 信道 状态 ， 换 言 之 ， 在 要 求 较 低 导 频 开销 的 情况 
SSA-CoSaMP 算法 才能 较 好 的 估计 出 
8(b) 可 以 得 ， 当 导 频 数 提高 到 750 时 ，SSA-CoSaMP 算法 略 优 


P.H 


言 道 状 态 。 从 图 


于 S-CoSaMP。 但 对 于 CoSaMP 算法 ， 导 频数 提高 到 750 仍然 


无 法 使 其 进行 有 效 运算 。 


图 8(o) 是 导 频 数 提高 到 850 的 仿真 
SSA-CoSaMP 算法 仍 优 于 S-CoSaMP $1 


图 。 当 SNR 低 于 15dB. 
E; 当 SNR 高 于 15dB 


E H- 2L E 86 H 
Ij, SSA-CoSaMP fil S-CoSaMP 算法 的 估计 性 能 基本 一 致 。 
101 
—Éé— S-CoSaMP ff i: 
一 一 SSA-CoSaMP 算 法 
一 6 一 CoSaMP 算 法 
10° 
wW 
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一 全 一 S-CoSaMP 算 法 
一 此 一 SSA-CoSaMP 算 法 
一 6 一 CoSaMP 算 法 


20 25 30 


(c) N, =800 


图 8 ju 


综合 本 节 的 仿真 数据 可 知 ，CoSaMP 算法 在 人 有 


算法 的 NMSE 性 能 对 比 


E 何 情况 下 性 


能 都 远 不 及 S-CoSaMP 算法 和 SSA-CoSaMP 算法 , 而 当 导 频数 
小 于 750，SSA-CoSaMP 算法 的 性 能 优 于 S-CoSaMP 算法 ， 并 
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H. SSA-CoSaMP WAI ERA BEME, 更 适用 于 稀疏 度 
未 知 的 实际 环境 。 实验 结果 表明 了 SSA-CoSaMP 算法 增强 了 稀 
玖 性 ， 同 等 估计 性 能 条 件 下 ， 降 低 了 导 频 开销 ， 节 省 了 频谱 资 
源 ， 在 低 导 频 开销 的 情况 下 ，SSA-CoSaMP 更 具有 优势 。 


6 ”结束 语 


本 章 根 据 大 规模 MIMO 信道 所 特有 的 空 时 共同 稀疏 性 的 
这 个 重要 特征 入 手 ， 从 动态 稀 下 度 自 适应 和 结构 化 两 方面 对 
CoSaMP 算法 进行 改进 ， 提 出 SSA-CoSaMP 算法 。 不 仅 优 化 了 
估计 性 能 ， 同 时 还 降低 导 频 开销 ， 节 省 能 源 消耗 和 频谱 资源 。 
符合 我 国 对 于 绿色 通信 成 为 发 展 趋势 的 要 求 。 仿 真 结 果 表明 ， 
无 论 是 低 导 频 , 还 是 低 SNR 的 状况 下 , 所 提出 的 方法 相对 传统 
基于 导 频 的 信道 估计 方法 都 有 明显 的 性 能 增益 。 

在 无 线 通信 环境 下 ， 结 构 化 特性 不 仅 存 在 于 实际 的 延迟 多 
径 域 中 ， 虚 拟 角 度 延迟 域 也 存在 稀 疏 度 结构 化 特性 的 特点 。 
此 , 接 下 来 的 研究 工作 主要 针对 大 规模 MIMO 天 线 阵 列 产生 的 
虚拟 角度 域 中 稀 玻 度 结构 化 的 问题 ， 使 结构 化 改进 方案 可 以 在 
虚拟 角度 域 中 应 用 ， 深 入 探究 结构 化 的 使 用 范围 ， 提 升 该 方案 
的 适用 性 。 
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